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Abstract 



We analytically and numerically investigate the 't Hooft-Bergknoff-Eller 
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mass of the lightest bound state. Its value is consistent with the precedent 
results. 
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I. INTRODUCTION 



Light-front (LF) quantization is believed to be an effective method for studying many- 
body relativistic field theories [0,0. The physical vacuum is equivalent to the bare vacuum 
in the LF coordinate, since all constituents must have non-negative longitudinal momenta 
defined by k + = (k° + k 3 ) / v2. This simple structure of the true vacuum enables us to avoid 
the serious problems which appeared in the Tamm-Dancoff (TD) approximation || in the 
equal time frame. Therefore, the TD approximation [|J is commonly used in the context of 
the LF quantization. 

The techniques have been developed [0-0] for solving LFTD equations in several models 
such as the massive Schwinger model ||, which is the extension of the simplest (1+1)- 
dimensional QED2 ||. Bergknoff []nj first applied LFTD approximation to the massive 
Schwinger model. He obtained the so-called Bergknoff equation, which is the light front 
Einstein-Schrodinger equation truncated to one fermion-antifermion pair. He obtained ex- 
cellent results for the lowest energy meson under 't Hooft's ansatz. He also discovered that 
it is necessary to include two fermion-antifermion pairs in order to study the excited states. 

Ma and Hiller || studied the lowest Bergknoff equation numerically. They devel- 
oped their numerical method for solving the Bergknoff equation using an idea based upon 
't Hooft's ansatz. 

Mo and Perry || suggested that even in the first excited state, most of the wavefunction 
consists of four fermion sectors when the fermion mass reaches zero. Therefore, four or more 
fermion sectors must be included so as to describe wavefunctions beyond the ground state 
bosons. 

In the bibliography |J , Mo and Perry presented an effective way to treat the ground state 
and the excited state in the massive Schwinger model. They concluded that to study the 
massive Schwinger model, the Jacobi polynomials are suitable as basis functions. Harada 
and his coworker [[?[] studied the massive Schwinger model with SU(2) flavor symmetry, in- 
cluding up to four fermion sectors. They used simpler basis functions, which are essentially 



equivalent to the Jacobi polynomials. Sugihara and collaborators |l I] numerically analyzed 
2-dimensional SU(Nc) Quantum ChromoDynamics(QCD) JT2], including four fermion sec- 
tors, by means of the basis functions of Harada et al. 

Although excellent papers exist concerning massless and massive Schwinger models and 
2-dimensional QCD ||i6|-|19|1 , including excited states, it is worth analyzing the '"t Hooft- 



Bergknoff-Eller" equation |U| , the extension of the 't Hooft-Bergknoff equation, in order to 
include both SU(N C ) and U(N C ) gauge theories. This is because there is a mathematical 
interest in the basis function method. There is no mathematical evidence that the con- 
ventional basis function expansion describes the wavefunction well; instead, the contrary is 
rather to be supposed, as there is evidence that the conventional method breaks down if we 
try to improve the approximation. We would therefore like to improve the basis functions 
so as to avoid such difficulties. 

Further, there is a so-called "2% discrepancy" problem, briefly summarized as follows: 
We expand dimensionless meson mass squared M 2 in terms of dimensionless quark mass m, 

M 2 = 1 + bim + b 2 m 2 + • • • , (1.1) 

where 
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M = — — and m 2 = -tt^. (1.2) 

Here, M is a mass of bound state, fh denotes a bare mass of quark, and n{Nc) stands for 
(N c 2 _ i _|_ a \ g i 

fi(Nc) = 77 • That is, we measure all masses in the unit of /x(iVc). Banks [13 

2%Nc 



et al obtained first order coefficient hi analytically, using the bosonization method, 

b x = 2e< E = 3.56214- ■•, (1.3) 



where 7^ = 0.57721 ■ • ■ is the Euller's constant. On the other hand, Bergknoff [nj found 
the value, 

61 = = 3.62759- •-, (1.4) 
v3 

which differs from Eq. (|1.3|) by 2%. In [P~5] , the authors suggested that the coefficient b\ = % 
was a variational invariant and that this discrepancy was ascribable to the contributions from 
the higher Fock sectors. Before proceeding on to consider the higher Fock sectors, we have to 
examine whether this discrepancy can be explained or not in terms of the lowest light-front 
Tamm-Dancoff equation using all possible basis functions. 

The 't Hooft-Bergknoff-Eller equation for two dimensional gauge theory is given in the 
form 

M 2 1>{x) = [\dyH(x,y)$(y) 

4(m 2 - 1) r\ $(v) fh 11 

where p denotes the finite part integral, a = 1 for U(Nc), and a = for SU(Nc). In Eq. (|1.5| ), 
we shifted the variable x total amount of — \ compared with the variable in [|l^|l^,[13|] , in 
order to show the symmetry of the wavefunction transparently. 

Mo and Perry concluded in that the Jacobi polynomials, (1 — 4x 2 )^P^ ,t3 (2x) in our 
notation, are the most suitable basis functions. This conclusion seems quite natural because 
the system of the Jacobi polynomials P^ ,l3 (2x) is an orthogonal complete set on the interval 
— I < x < I with respect to the weight function (1 — 4x 2 ) /3 . Harada and collaborators 
suggested using the simpler basis functions, (1 — 4x 2 )^ +J and x(l — Ax 2 ) l3+ ^ in our notation, 
which are equivalent to the abovementioned ones. According to Harada et al, one can expect 
that the wavefunction could be expanded as follows: 



<3>(x) 



N 



lim Y aAl-Ax 2 f +j 

N 

lim Y hxtt -Ax 2 f +j . 

3=0 



(1.6) 



Here, we have used the fact, as is shown in Appendix that the Eller equation does not 
mix even and odd functions with each other. 
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The exponent (3 and the quark mass m are related to each other by the equation |T2|JT0[| 

(m 2 -l)+/3vrcot/37r = 0. (1.7) 

The authors of references |||l(],[7|,[ll| adopted the positive smallest solution /5 (m) of Eq.( |1.7|) 
as (3 in Eq.( |1.6| ). That is, for small m, 

Mm) = ^(i-?l) + 0(m 5 ) = A/m + A) 3 ™ 3 + ■■■. (1.8) 
7r y 10 J 

In the following sections, we try to determine the coefficients a n 's, according to the prede- 
cessors. 



II. CONVENTIONAL BASIS FUNCTION METHOD 



In this section, we restrict ourselves to the case a = where the gauge group is SU(N C ). 
Mo and Perry, and Harada and his collaborators, presented effective way to determine the 
coefficients. We will briefly reproduce their procedures. By the use of the expansion in 
Eq.(O) truncated to given finite number N for the wavefunction $, we multiply both sides 
of Eq. (|1.5| ) by (1 — 4x 2 ) /3 ~ H and integrate them over x, then we obtain 

(2.1) 



M 2 Na = Ha, a = *[ao, ai, • • • , a n -i\- 
Here N and H are n x n matrices and are given by 



i dx(l — Ax 



and 



Hi 



4(m 2 



2k 1 ' 2 (m 



+ 



\ dx(l - Ax 2 ) 2 ^- 1 

T{2f3 + i+j) 
>T{2f3 + i + 3 + l/2) 



P 



2T(2f3 + 3/2) ' 

1 -4x 2 )P +i (l -4y 2 )P + i 



(2.2) 



dxdy- 



(y-x) 



2(3 + i+ j 



B(j3 + i,/3 + i)B(J3 + j,p + j). 



(2.3) 



See Appendix of 0. So-called "norm" matrix N appeared in the above equation because the 
basis functions we have used are not orthonormalized. In order to have eigenvalues of the 
generalized eigenvalue equation given in Eq. ( |2.1| ), we have to solve the eigenvalue problem 
for norm iV first, i.e., 



Nm = XM. 



Next, we introduce a transformation matrix W by 



W 



A7' 



(2.4) 



(2.5) 
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Then, we can transform Eq. (|2.1|) into a usual eigenvalue problem of the form 



M 2 b = l WHWb, a = Wb. (2.6) 

We can solve Eq. (|2.6| ), numerically. For N = 3 and m = 0.01, we find, for the ground 
state boson, 

P = 0.00552328, M 2 = 0.0366342, 

a = 1, ai = 0.00203562, a 2 = -0.000579369, a 3 = 0.000165813. (2.7) 

The values of the LHS and the RHS of Eq.( |1.5| ) are shown in Fig. 1. The coincidence of the 
LHS and the RHS is high for small values of x. For x ~ ±1/2, the behavior of the LHS and 
the RHS are quite different. There are sharp spikes at the end points. This behavior is not 
changed much even if we improve the order of approximation. 

Note here that in order to solve the generalized eigenvalue problem, the norm matrix 
should be positive definite. We cannot advance the above procedure beyond iV ~ 12, 
because some of the eigenvalues of the norm matrix N become almost zero or negative. We 
will examine mathematically this approximated wavefunction in detail in the next section. 



Fig. 1 



III. AN INSPECTION OF THE CONVENTIONAL BASIS FUNCTION METHOD 

A. Behavior of wavefunction around x = 

We introduce linear map C by 

C: f ^ Cf such that (Cf) (x) = p [* dy ^ X ^ V \ (3.1) 

j-\ y 2 

After some tedious but not particularly difficult calculations, we find that 



C : (1 - 4*y ~ -^ 1/2 r(p + l/2) F ^ 1/2 ~ ^ 1/2; A%2) 
C:x (l- 4*Y „ ~ 8 * 1/2 + X F(2 ' 1/2 ~ 3/25 



(3.2) 



where F(a,b;c;x) is the Gauss' hypergeometric function [ S0[ . See Appendix A. 

We restrict ourselves to the case where the wavefunction is an even function, because we 
are interested in only the ground state meson. Thus, we are led to 

N 

lim Y a n M 2 (l-Ax 2 f +n 
7V - +O °^o 



lim a 

n=0 



r 4(m 2 - 1)(1 - 4x 2 )^ + ** 1/2 ry + Xim F{1 ' V2 " P " n; 4X2) 



. (3.3) 
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Now, we can examine whether the numerical result obtained so far satisfies Eq. (|1.5|) or 
not. We substitute Eq. fl2.7Q into the above equation, and obtain 

LHS(x) = 0.0366936 - 0.00101207a; 2 - 0.00165597 x 4 - 0.00466671 x 6 + 0(x 8 ) (3.4) 

and 

RHS(x) = 0.0367133 - 0.0038283 x 2 + 0.0503141 x A - 0.210419 x 6 + 0{x 8 ). (3.5) 

Thus, the LHS and the RHS coincide, within numerical errors, with each other only up to 
O(x ). We calculated coefficients a n 's up to n = 10, but the coincidence between the LHS 
and the RHS is not much improved. See Fig. 1. 



B. Behavior of wavefunction near x 



±1 



±| is 



Bergknoff suggested that the behavior of wavefunction $ near end points x 
important in order to calculate the mass eigenvalue M 2 . In fact, Eq.( |1.7| ) is derived by 
demanding that the most singular part, that is, the coefficients of (1 — 4x 2 ) /3 ~ 1 on the 
RHS in Eq. fll.5p must be cancelled, as there is no such term on the LHS. According to 
Bergknoff's suggestion, we will examine the behavior of wavefunction near the end-points 



beyond O (Jl — Ax 2 ^ j . We set 4x 2 = 1 — e. By the use of the identity for the hyperge- 

ometric functions, that is, Eq.(15.3.6) in [20], which is valid for a + b — c ^ integer, we can 
expand the RHS of Eq.(EO) around e = and have 



LHS = M 2 Y, a n e p+n . 

n=0 



(3.6) 



and 



RHS = 4(m 2 - 1) ]T a n e p+n - 1 

n=0 

oo n 

+ 47r£5> i (/3 + j)cot (tt (p + j)) 



•1/2), 



n=0 j=0 



2n^ 2 e 



n=0 



3=0 



(n-j)\ 

T(P + j + l)(l/2-(3-j) n 
(/3 + i - l)r(/3 + i + 1/2) (2 - /3 - i) T 



(3.7) 



Substituting the numerical solution, which is given in Eq. fl2~T|) , into Eqs. (|3.6|) and (p.7|), 

we have 



LHS(e) = 0.0366342 + 7.45734 x 10~ 5 e^ 1 - 2.12247 x 10~ 5 e p+2 + 0(e 



and 



RHS{e) = 0.551656 - 0.525734 e p + 2.09072 e - 2.08014 e^ 1 + 0(e 2 ). 



(3.8) 



(3.9) 



Only e 13 1 terms on the LHS and the RHS coincide with each other, because we define (3 
so that the coefficients of e^ -1 cancel each other on the RHS. Note that the LHS does not 
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contain e n terms with non-negative integer n, while the RHS does. These tendencies are not 
changed even if we calculate a n 's for n=3, 6 and 10. 
If we rewrite the first two terms of Eq. ( |3.9| ) as 

(3.10) 



0.551656(1 - e 13 ) + 0.025922e /3 , 



we can see the origin of the spikes at the endpoints in Fig. 1. The spikes arise from the 
existence of the constant term in the wavefunction. The wavefunction $(x) = 1 is the exact 

and M 2 = in SU(Nc). Thus, one may expect that the 



m 



solution of Eq. (|1.5|) for 

spikes at the end points are closely related to the existence of the massless bound state in 
SU(Nc). This is not the case, however, because we can easily see that the constant term, 
in the wavefunction, is allowed if and only if m 2 = 0. We may conclude that the spikes are 
nothing but the artifact which arose from the fact that we have used the improper basis 
functions in Eq.( |1.6| ). In fact, we may remove the spikes if we use the suitable wavefunction. 
Refer to the next section. 

The basis function given in Eq.( |1.6| ) cannot be a good mathematical approximation of 



the true wavefunction. The reason is as follows: If we truncate the series Eq. ( |1.6|) to N, we 
expect that Eq.fll.5p holds up to 0( 



We have only N + 1 parameters ai, a 2 , 



and M 2 . On the other hand, we have 2N equations up to 0(e /3+7V_1 ). No consistent solution 
can exist in general in this case. 

The main difficulty comes from the fact that C maps (1 — Ax 2 ) 13 not only to the terms 
(1 — 4x 2 ) /3+ ' 7 ~ 1 with non-negative integer j but also to the terms (1 — Ax 2 )i . One may expect 
that the above difficulty is avoidable if the additional terms (1 — 4x 2 ) J are introduced in 
Eq.( |1.6| ). However, we can easily see that coefficients of such terms must cancel. By the use 
of the identity, that is Eq. (15.3.11) in |20| 



we have for n > 2 



C:(l-4x 



2\n 



-L J i, n 



2(n-l)ft (2-n) k 

ttI/2 oo ( n ) (-1/2) 

• 2\n-l \ - V Jk \ 1 Jk n 



[-l + Ax 



E 



r(l/2-n) v '" ' fo k] ( k + n - !) ! 

log(l - Ax 2 ) - ip(k + 1) + ip{-l/2 + k) 



(1 -Ax 2 f x 



(3.11) 



Here, ip{z) denotes the digamma function. An analogous formula for n = 1 holds. See 
Eq.(15.3.10) in [^DJ. If we introduce a term (1 — Ax 2 ) n in Eq. fll.6p with positive integer n, a 
new singular term like e n ~ l loge appears only on the RHS. Note that the exponent n of e is 
the same as that of the introduced term. Thus all the coefficients of (1 — Ax 2 ) n should be 
zero. 



IV. NEW BASIS FUNCTION 

We must notice that there are infinite solutions of Eq. fll . 7| ) in addition to the solution 
jiven by Eq. (|1.8|). In fact, we see that 
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/3 n (m) = n + 1/2 - - ■ — - - — a + O 



+ m 2 



(n + l/2)7r 2 3(n + l/2) 3 vr 4 \(n + l/2) s ir* 



+ 0(m 4 



> + 1/2)tt 2 3(n+l/2)V V( n + W*)*** , 
= p n ° + p n 2 m 2 + ---, n = 1,2,3 (4.1) 



From Eqs. fll.8Q and (|4.1|) , we are led to 

< /3„(m) - /3 (m) - n < 1/2, 

< p n {m) - p k {m) - {n - k) < 1. (4.2) 

The above relations imply that Eq. (|1.6| ) never incorporates terms like (1 — 4x 2 )^ n+J with 
positive integer n and non- negative integer j. 

We posit that the wavefunction is given by an infinite series 

N N-n „ I \ , • 

= Jim^E E C « J (l - 4x 2 )^ (m) ' . (4.3) 

n=0 j=0 

For counting the number of free parameters and the number of nontrivial equations, we 
consider the truncated wavefunction to given finite N. The truncated wavefunction includes 
the term (1 — 4x 2 ) /3jv and all the other lower order terms. We require that Eq.( |1.5| ) should 
hold up to 0(1 — 4x 2 ) /3jv ~ 1 . For each given value of m, the unknown parameters are M 2 and 
cj except for cq° = 1. Thus the number of the parameters is ( Af + 1 K jV + 2 ) _ On the other hand 
we have jV( -^ +3 ' > nontrivial equations. See Table |I[ 



Table 1. 



The number of parameters is larger than that of non-trivial equations by 1. Thus, we 
can solve the equations for cj in terms of M 2 . Another equation of use to us is obtained 
by multiplying both sides of Eq.( |1.5| ) by $(x) and integrating them over x, 



i ii 
M 2 J \ dx\${x)\ 2 = J\ j\dxdy${x)H{x,y)${y). (4.4) 



It should be noted here that Eq.( |4.3|) is, mathematically, the most general expansion. 
This means that there is no room to introduce any other additional terms like d(l — 4x 2 ) 7 
for 7 7^ /3 n + j with non- negative integers n and j. If we introduce such terms, then the 
following equality should hold 

= 4rf (m 2 - 1 + 7T7 cot(7T7)) (1 - 4a; 2 ) 7 " 1 . (4.5) 

This demands that d = 0. 

In the following subsections, we will examine our new basis function in detail both 
analytically and numerically. 
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A. An analytic approach 



In this subsection, we will restrict ourselves to the case where N — 1. Up to 0{e^ l ) 
we have four equations: 



gft-i . o = 4( m 2 _ i) Cq o + 47T/?oCot(7r/5 )c , 

i r(2 + /5 ) 



(4.6a) 



(i-A))r(i/2 + A)) Co /3 r(3/2 + /3 ) 
r(i + ft) 

(l-^)r(l/2 + y9i) 

' o r(i + A,) rl r(2 + /3 ) (l + A) 



+Ci 



a 
+ 4 



r(3/2 + /3 ) u r(| + /5 ) A r(3/2 + ^) 



: c °M 2 = 4c 1 (m 2 - 1) + 4vr cot(vr/? ) 
e^- 1 : o = 4(m 2 - l)ci° + 4tt/3i cot( 7 r/3 1 )ci Q . 



c °/3o 



+ (l + /3o)c 1 



(4.6b) 

(4.6c) 
(4.6d) 



Equations (|4.6a|) and (|4.6d|) are automatically satisfied. Since c ° = 1, we can solve 
Eqs. (|4.6b|) and ( |4.6c| ) for c Q l and ci° in terms of m, M, @ , and (3\. 

Now, in order to solve the above equations for M, we assume that all physical quantities 
can be expanded in terms of quark mass m. That is, 

M 2 = b + hm + ■ ■ ■ . (4.7) 

Thus, we have, up to 0(m), 

i 

dx\$(x)\ 2 =< $|$ > 

2 

2/V (22 + 8/3x (1 - 3 log(2)) - 12 log(2) + a (2 - 6 log(2) - px° (5 - 6 log(2)) 



1 + m 



x 



3 (-2 + a (-1 + A ) - 4ft ) 



&im 



2 (l - A°) 



2 + a + 4ft - aft ( 



(4.8) 



and 



2 / 2 



/ i dxdy§(x)H(x,y)$(y) =< $|#|$ 



> 



a + m 



2 (l + 2/3i°) (3 + /V V) + a (3 (-1 + ft ) + (ft 1 )' (48 + 6ft - tt 2 



+ ftV - 361og(2) - 36ft°log(2) 
2 (l - ft°) 



3ft 1 (-2 + a(-l+ft°)-4ft ) 



+ 61m 



2 + a + 4ft u - a/?/ 



(4.9) 



We may expect that coefficient 61 depends on a, as both coefficients of m in Eqs.([4.8|) 
and ( |4.9|) explicitly depend on a. However, this is not the case. Indeed, from Eq. Q4.4|) , we 
are led to 
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M 2 = a + I + m + °(™ 2 ) 




(4.10) 

Note that we obtained the first line in the above expansion without referring to the explicit 
value of Pi. We did not reproduce the result of Banks et al, but that of Bergknoff. The ap- 
proximated wavefunction in case m = 0.01 and iV = 1 is shown in Fig. 2. The approximation 
used here is so rough that the coincidence of the LHS and the RHS is poor. Nevertheless, 
the behavior of the RHS near the end points is quite calm compared with the results of 
the conventional basis function method. The smoothness of the RHS is quite natural. As 
mentioned previously, the wavefunction given in Eq. ( [4.3| ) is the most general. If we truncate 
the wavefunction up to order N, it becomes smooth. So we may also expect the RHS of the 
't Hooft-Bergknoff-Eller equation to be smooth. 



Fig. 2 



We may expect that the coincidence of the LHS and the RHS will be improved if the 
higher order terms are included. In cases where iV > 2, we cannot treat things analytically. 
We will attempt to solve Eq. fll.5p numerically by the use of new basis function in the next 
subsection. 



B. A Numerical approach 



In general, we have 



= MHn.r) - r x dyH(x,y)$(y) 

"2 

OO 



l-4x 1 =e 



-4 2 c n ° (m 2 - 1 + np n cot irp n ) i 

n=0 

oo oo / j+1 (- 1 / 2 ) 

- 53 53 M V - 4(m 2 - l)c^ +1 - 4tt 53 c n k (P n + k) cot npj ^ 

n=0 j=o y 



k=0 



(j + l-k)\ 



J3n+j 



+ 



7T 



1/2 oo oo 



AT(P n +j + l) 



aV(p n + j + l) 



* n =0j=0 

oo / oo oo 

-S^E EE C " 
k=0 \n=0j=0 



(P n + j - l)T(P n + j + 1/2) T(P n + j + 3/2) J 

T(P n+J + l) { l / 2 -^-j) k+ i 



(4.11) 



(P n+j -l)T(p n +j + l/2) (2-P n -j) 

\ / fe+i 

Here, e = 1 — Ax 2 as before. Of course, the first line in Eq.( (4.11|) cancels automatically 



because of the definition of /3 n 's. Then, suppose that we truncate series in Eq.( |4.3|) to 
0(e^ N ). That is, we set cj = for n + j > N. For a given m, we put M 2 = M 2 . 
We can then solve Eq. (|4.11| ) for cj in terms of M,. We thus obtain the Mj dependent 
truncated wavefunction, say, $(x; Mi). We can calculate a new mass eigenvalue M i+ i using 
this wavefunction as 
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<HM i )\H\HM i )> 
"' _J " < $(Mi)|$(Mi) > ' 



We can use Eq.( [4.10D as Mq. For iV < 15 and m = 0.01, mass M converges in 5 iterations 



For < m < 0.5, we obtain M 2 's which are summarized in Table |TJ. We can fit them by 
polynomials, as follows: 

M 2 (a = 0, m) = 3.62763m + 3.58027m 2 + 0.0683573m 3 + 0(m 4 ) 

M 2 (a = 1, m) = 1 + 3.62421m + 3.34492m 2 + 0.213839m 3 + 0(m 4 ). (4.13) 

It should be noted here that the coefficients of m are consistent with Eq. fl4.10Q and 
Bergknoff's result. In order to see the efficacy of this new basis function expansion, we 
sho w the wa vefunctions in Fig. 3. 

Table 2 

Fig. 3 I 



V. SUMMARY AND DISCUSSION 

In the preceding sections we have introduced the new basis function and calculated the 
mass eigenvalue of the bound state using the new basis function. We have found that (1) 
the new basis function gives an effective approximation of the wavefunction, and (2) the 
mass eigenvalues are consistent with the results of the precursors. In the remainder of this 
section, we will discuss the 2% discrepancy problem. 

Let us consider the wavefunction given by 

$(x) = (1 - 4x 2 ) 70 + ai(l - 4x 2 ) 71 , (5.1) 

where we assume only 70 = 7o 1 m + 7i 3 m 3 + 0(m 5 ), 71 = 71 + 71/m + 0(m 2 ), and a\ = 
a\ l m l / 2+s + O(m). We are then led to 

M 2 = < $|if 1$ > / < $1$ > = a + [Xr + m + 0(m 1+2S ). (5.2) 

\V 3 / 

This relation holds independently of the details of 70, 71, and ai, except for certain assump- 
tions which were made before Eq. (|5.2j ). The coefficient of m in Eq. (|5.2|) has the minimum 
value ^| when 70 1 = Po 1 = We may therefore conclude that Eq. ( |5.2| ) holds universally, 
provided that: 

(1) The wavefunction can be expanded as a power series of (1 — 4x 2 ), like 

= (1 - 4X 2 ) 70 + ET=i - 4 ^ 2 ) 7j , 7o < 7i < • • ■, 

and 

(2) The coefficients of the series, Oj's, are of order m 1 / 2+<5 . 



An almost identical result has been obtained by Harada et al [I5|, in which they have 
restricted themselves to a case where 70 = /3 , and jj = /3 + j. Our conclusion is a 
generalization of Harada et aPs result. We, therefore, cannot solve the "2% discrepancy" 
problem in the context of the 't Hooft-Bergknoff-Eller equations. 
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APPENDIX A: PROOF OF EQ.(3) 



We will prove Eq. fl3.2|) . For monomial x n , by the use of the definition of the finite part 
integral, we obtain 



Here. 



Ax n 1 — 2x 

C:Xn ^ ~l^Ax~ 2 + UX bg T+2x~ + UX) 



(Al) 



m^k I / 1 



k=2 
( n-1 

E 



X 

2 



fc-i 



1 

x 

2 



k-V 



x 



m—k 



{2k + l)x 



2k 



k= - {2(n - k) - 1}2 
" _J (2k + 2)x 2k+1 



2n-2k-2 ' 



lS{2(n-fc)-l}2 



2n-2fc-2 ' 



for m = 2n, 



for m — 2n + 1. 



(A2) 



Using identities 



Ax 



2n 



n=0 



Ax 2 



in— fc+l 2n 



n=0 fc=0 



(A3) 



1 - 2x 



Y 2n c 2n x 2n ~ l log = -Y 

r-i l + 2x f-i 



n=0 



n=0 



.£ C2fc 2(^-fc) + l 



and 



E °2nf 2n(x) = ~Y 



n=l 



n=0 



™ (2n + i)2 2 "" 2fc+1 
^ 2(n-fc) + l 



2d 



we see that, for a given even function Y^^Lo c 2n,x 2n , 



oo An—k+1 



n=0 n=0 fc=0 Z V n — K ) + 1 



Power series expansion 



:i-4*y = £ 



-/3)„(4z 



2\n 



n=0 



(A4) 



(A5) 



(A6) 



(A7) 
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Y(a + n) (—B) n A n 

holds where (a) n = — — is the Pochhammer symbol. Thus, for c 2 „ = p — 

r(a) ra! 

C2n+i = 0, we see that 

^ : (1 - 4*Y - - 4 r(i/ 2 ( ^f ^(^ V2 - # 1/2; 4* 2 ), 
where F(a, 6; c; z) is the Gauss' hypergeometric function. Analogously, we have 

C:x(l- 4 X y - J t(l/2 + P) xF ^ 1/2 ~ # 3/2; 4x2) - 



13 



REFERENCES 



[1] S. J. Brodsky, H-C. Pauli, and S. Pinsky, Phys. Rep. 301, 299 (1998). 
[2] M. Burkardt, Adv. Nucl. Phys. 23 1 (1996). 

[3] I. Tamm, J. Phys. (Moskow) 9, 449 (1945); S. M. Dancoff, Phys. Rev. 78, 382 (1950). 
[4] R. J. Perry, A. Harindranath, and K. G. Wilson, Phys. Rev. Lett. 65, 2959 (1990). 
[5] Y. Ma and J. R. Hiller, J. Comp. Phys. 82, 229 (1989). 
[6] Y. Mo and R. J. Perry, J. Comp. Phys. 108, 159 (1993). 

[7] K. Harada, T. Sugihara M. Taniguchi and M. Yahiro Phys. Rev. D49, 4226 (1994). 
[8] S. Coleman, R. Jackiw, and L. Susskind, Ann. Phys. (N.Y.) 93, 267 (1975); S. Coleman, 

Ann. Phys. (N.Y.), 101, 239 (1976). 
[9] J. Schwinger, Phys. Rev. 128, 2425 (1962). 
[10] H. Bergknoff, Nucl. Phys. B122, 215 (1977). 

[11] T. Sugihara, M. Matsuzaki and M. Yahiro, Phys. Rev. D50, 5274 (1994). 
[12] G. 't Hooft, Nucl. Phys. B75, 461 (1974). 

[13] K. Hornbostel, S. J. Brodsky and H.-C. Pauli Phys. Rev. D41, 3814 (1990). 
[14] T. Banks, J. Kogut, and L. Susskind, Phys. Rev. D13, 1043 (1976). 
[15] K. Harada, T. Heinzl, and C. Stern Phys. Rev. D57, 2460 (1998). 
[16] S. Huang, J.W. Negele and J. Polonyi, Nucl. Phys. B307, 669 (1988). 
[17] B. van de Sande, Phys. Rev. D54, 6437 (1996). 

[18] A. Duncan, S.Pernice, and E. Schnapka, Phys. Rev. D55, 2422 (1997). 
[19] O. Abe, G. J. Aubrecht, and K. Tanaka, Phys. Rev. D56, 2242 (1997). 
[20] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, (Dover Publ., Inc., New York, 1972). 



14 



TABLES 



TABLE I. Number of nontrivial equations 



terms 


range of n 


range of j 


number of nontrivial equations 


(1 - 4x 2 ) /3 "" 1 


~ N 




(automatically satisfied^ 


(1 - 4xY n+i 


~ N 


~ N- n- 1 


N(N+1) 
2 


(1 - 4x 2 y 




0~ N- 1 


TV 


TABLE II. Numerical results for bound state mass M 2 
function of quark mass m 


in SU(Nc) and U(N C ) models as a 


m 


0.01 


0.10 0.20 


0.30 0.40 0.50 


M 2 in SU(N C ) 


0.036634 


0.398634 0.869282 


1.412358 2.028271 2.717415 


M 2 in U(N C ) 


1.036607 


1.396177 1.860377 


2.394130 2.998685 3.675073 
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FIGURES 




FIG. 1. The comparison of the relative values of both sides of Eq,([I.5|) for a = 0. The wavefunc- 
tion in Eq,( |1.5| ) was approximated by Eq.(|l,6|) with N = 3 and Eq,( |2.7| ). The solid line represents 
the LHS and the dotted line stands for the RHS. The RHS with N = 9, which is indicated by the 
dashed line, is exhibited for comparison. 



1 -05| — . 1 1 , , 1 , , 1 , r 




FIG. 2. The comparison of the relative values of both sides of Eq. (|l.5| ) for a = 0. The wave- 
function in Eq.([1.5|) was approximated by Eg. ( [4. 3D with N = 1 and Eq.( |4.1C| ) with m = 0.01. The 
solid line represents the LHS and the dashed line stands for the RHS. 
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1.02, — . , 1 , , 1 1 , 1 1 r 




X 



FIG. 3. The convergence of the new basis function expansion for m = 0.01 fixed. The thin solid 



line represents the LHS in Eq.(L5), provided that the wavefunction was approximated by Eq.( [4.3[) 
with N = 15. The dotted line denotes the RHS with N = 2, the dot-dashed line exhibits the RHS 
with N = 3, the dot-dot-dashed line represents the RHS with N = 4, the dot-dash-dashed line 
stands for the RHS with N = 5, and the dashed line exhibits the RHS with N = 10. The thick 
solid line indicates the RHS in Eq. (|1.5| ) with wavefunction given in Eq.([4.3j) with N = 15. 
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